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ABSTRACT 

The morphological properties of galaxies between 21 mag < / < 25 mag in the Hubble 
Deep Field are investigated using a quantitative classification system based on mea- 
surements of the central concentration and asymmetry of galaxian light. The class 
distribution of objects in the Hubble Deep Field is strongly skewed towards highly 
asymmetric objects, relative to distributions from both the HST Medium Deep Survey 
at / < 22 mag and an artificially redshifted sample of local galaxies. The steeply rising 
number count-magnitude relation for irregular/peculiar/merging systems at / < 22 mag 
reported in Glazebrook et al. (1995a) continues to at least / = 25 mag. Although these 
peculiar systems are predominantly blue at optical wavelengths, a significant fraction 
also exhibit red U — B colours, which may indicate they are at high redshift. Beyond 
Glazebrook et al.^s magnitude limit the spiral counts appear to rise more steeply than 
high-normalization no-evolution predictions, whereas those of elliptical/SO galaxies only 
slightly exceed such predictions and may turn-over beyond I ^ 24 mag. These results 
are compared with those from previous investigations of faint galaxy morphology with 
HST and the possible implications are briefly discussed. The large fraction of pecu- 
liar/irregular/mcrging systems in the Hubble Deep Field suggests that by / ~ 25 mag 
the conventional Hubble system no longer provides an adequate description of the mor- 
phological characteristics of a high fraction of field galaxies. 
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1 INTRODUCTION 

The Hubble Deep Field (HDF) is a four square arcminute 
area of sky at RA 12h 36m 49.4s Dec +62° 12' 58.0" (J2000) 
observed during December 1995 with the Hubble Space Tele- 
scope (HST). The images, obtained during 150 orbits with 
the Wide Field and Planetary Camera 2 (WFPC-2), rep- 
resent the deepest optical imaging survey yet undertaken 
(Williams et al. 1996). Four broadband filters were used 
whose wavelength coverage samples the range between the 
ultraviolet and near infrared (roughly UBVT). In the I- 
band (F814W), the HDF data reaches ~ 3 mag fainter than 
the deepest ground-based observations at similar wavelength 
(Small et al. 1995), and ~ 1 mag fainter than the deepest 
7-band observations undertaken previously with the HST. 
The position of the HDF was pre-selected to be a "typical" 
sample of the high galactic latitude sky in terms of galaxy 
surface density, while also avoiding bright stars or other ob- 
jects whose presence might compromise the depth of the 
observations. 

A major objective of the HDF program is to determine 



the nature of the rapidly evolving faint blue galaxy popu- 
lation seen in ground-based spectroscopic and photometric 
surveys conducted over the past decade (Broadhurst, Ellis, 
& Shanks 1988; Tyson 1988; Colless et al. 1990; Ellis 1990; 
Koo & Kron 1992; Lilly 1993; Tresse et al. 1993; Cowie et al. 
1994; Glazebrook et al. 1995b). Many of these objects are 
distant {z > 0.3), and currently only the HST provides the 
resolution needed to characterise their morphological prop- 
erties. Results to date have proved difficult to reconcile how- 
ever, in part perhaps because of the inherently qualitative 
nature of galaxy classification. Number counts from a sam- 
ple of 301 galaxies with / < 22 mag taken horn 13 WFPC-2 
fields from the Medium Deep Survey (MDS) have been re- 
ported by Glazebrook et al. (1995a). These authors classified 
galaxies by eye and found the number-magnitude relation- 
ships for elliptical and spiral populations to be consistent 
with the predictions of high-normalization no-evolution cos- 
mological models, but that the number of irregular/peculiar 
(possibly merging) galaxies were at least an order of mag- 
nitude greater than predicted by the same models. A sim- 
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Hubble Deep Field (21 < 1 < 25) 
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Figure 1. Asymmetry vs. central concentration for galaxies in the Hubble Deep Field. RSE's visual morphological classifications are 
keyed to the plot symbols; E/SO's are shown as ellipses, spirals earlier than Sd are shown as spirals, and irregulars/pcculiars/mcrgcrs are 
shown as asterisks. The sectors subdivide the diagram into regions where each of these morphological types dominates. Representative 
error bars are shown. The dotted polygon shows the "convex hull" enclosing the artificially redshifted local sample of objects described 
in the text. 



ilar conclusion was reached by Driver et al. (1995) based 
on number counts from a single deeper MDS WFPC-2 field 
to 7=24.25 mag and by Abraham et al. (1996) from auto- 
mated morphological classification of 521 galaxies from a 
larger number (21) of WFPC-2 MDS fields to 1=22 mag. 

In contrast, HST observations reported by Schade et al. 
(1995) of a small but spectroscopically complete sample of 32 
galaxies from the Canada-France Redshift Survey (CFRS) 
indicate that the morphological mix at z ~ 0.75 is similar to 
that seen locally, except for a population of "blue nucleated 
galaxies" (BNGs) which Schade et al. suggest comprise ~ 
35% of the faint blue galaxy population. The existence of 
another major new population of morphologically distinct 
objects, linear "chain galaxies" that may be proto-galactic 
systems, has recently been reported by Cowic et al. (1995) on 
the basis of HST observations of faint objects (7 > 23 mag) 
identified in the Hawaii Deep Survey. 

A major advantage of the HDF dataset over the ear- 
lier images of Glazebrook et al., Schade et al., Driver et al., 
and Cowie et al. is the improved signal-to-noise level aris- 
ing from the considerably longer HDF exposure time. The 
exposure times available to those authors were 5-8 hours 
in F814W, which is substantially less than the ~35 hours 
devoted to the HDF in this filter. The improved signal to 
noise at faint limits enables us to extend the work already 
begun with the MDS, whose wider spatial coverage nicely 
complements the HDF data. 

Although interesting questions (which the HDF may ul- 
timately answer) concern the very faint (7 >25 mag) galaxy 



counts, in this short note we are concerned with understand- 
ing the morphology and colour distributions of the ~ 300 ob- 
jects brighter than 7=25 mag. As described below, at this 
limit the bulk of the galaxy population is detectable in all 
four HDF passbands, and morphological classifications are 
as robust as was the case in our earlier analyses of the MDS 
data to 7=22 mag. 

A plan of the paper follows. In §2 we review our mor- 
phological classification techniques in the context of the 
HDF data. Section 3 discusses the colour distribution as a 
function of morphology and, in §4, we combine our dataset 
with that of the MDS, enabling us to extend Glazebrook 
et al. ( 1995a) 's morphological galaxy count diagram three 
magnitudes fainter. A number of new trends appear whose 
implications are briefly discussed. 



2 MORPHOLOGICAL ANALYSIS OF THE 
HUBBLE DEEP FIELD 

Working with the "dithered" (rcsamplcd by a factor of 2.5 
and stacked) version 1 HDF images released by STScl, 
we constructed a source catalog using the APM faint- 
object photometry software (Irwin 1985). At very faint 
(7 ^ 27 mag) levels, the number of objects in the result- 
ing catalog depends sensitively on the smoothing scale, de- 
tection threshold, and minimum object area used by the 
softwaxe, but in the regime where morphological classificar 
tion is reliable (7 < 25 mag) the final catalog is robust. The 
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photometry was calibrated on the instrumental magnitude 
system defined by Holtzmann et al. (1995). Glazebrook et 
aZ.(1995a) used the same 7-band calibration in their MDS 
analysis and thus the results are directly comparable. The 
structural properties of galaxies on the J-band frame in the 
magnitude range 21 mag < 7 < 25 mag were investigated 
with the automated morphological analysis system used by 
Abraham et al. (1996) to characterise objects from the MDS. 
Visual classifications were also made by two of us (RSE and 
vdB) using the classification system discussed by Glazebrook 
et al. (1995a). The reader is referred to these earlier papers, 
and to Abraham et al. (1994), for further details of these 
procedures beyond the outline given hero. 

In the automated classification procedure, the galaxy 
images were first isolated from the sky background by ex- 
tracting contiguous pixels at least la above the sky level. 
For each galaxy an asymmetry index. A, was meeisured by 
rotating and self-subtracting the galaxy imago from itself. 
A second parameter in the classification procedure, central 
concentration C, is determined from measurements of the 
intensity- weighted second order moments of the galaxy im- 
ages. Asymmetric or disturbed galaxies generally have high 
values of A and low values of C, ordinary spirals have in- 
termediate values of both parameters, and early-types are 
distinguished by high central concentration. 

In order to compare directly the measurements of C 
and A obtained from "dithered" HDF images to those ob- 
tained from the MDS, the 23 galaxies in the HDF with 
21 mag < / < 22 mag (the magnitude limit of the Medium 
Deep Survey) were rebinned to the original sampling of the 
WFPC-2 camera and degraded by adding Poisson noise 
to simulate the appearance of these objects under typi- 
cal 5000s) MDS exposures. A comparison of measure- 
ments from the original and degraded images indicated that 
the HDF dithering procedure introduces a small system- 
atic change in the measured asymmetry value of the HDF 
data, = 0.053 ± 0.02, relative to undithered MDS data. 
This small ofi^set was therefore subtracted from the measured 
asymmetry values when comparing the HDF and MDS. It 
is emphasized that applying this correction decreases the 
the total number of asymmetric objects in the HDF. Ran- 
dom errors were determined by Monte Carlo simulation (c.f. 
Abraham et al. 1994 and 1996) and are no larger than 
AC = AA = 0.07 to 7=25 mag. Both C and A depend 
on the rest frame isophote level at which the measurements 
are made. The la limiting isophote of the HDF data (~ 26.5 
mag arcsec"^) is 1.5 - 2 mag arcsec"^ deeper than the cor- 
responding limiting isophote for more typical HST data ( eg. 
from the MDS), so bright galaxies in the HDF are probed to 
larger radii than corresponding objects in the MDS. Strong 
(1 -|- z)'^ cosmological dimming is expected to roughly syn- 
chronize the rest-frame limiting isophotes for the MDS and 
HDF samples at faint magnitudes, but a detailed assessment 
of this effect must await redshift information (discussed be- 
low). 

The distribution of HDF galaxies in the log (A) vs. 
log(C) plane is shown in Figure 1. The visual classifications 
of RSE are indicated by the plot symbols on this Figure, and 
are in excellent agroomont with classifications based on posi- 
tion on the log{A) vs. log(C) diagram. Using the visual clas- 
sifications and the simulations described below, three sectors 
were defined on the diagram, subdividing the galaxy pop- 



ulation into irregular/peculiar/mcrgors, spirals, and E/SO 
systems (the same bins adopted by Glazebrook et al. 1995a, 
and Abraham et al. 1996). The HDF asymmetry values are 
skewed markedly towards highly asymmetric galaxies rela- 
tive to the MDS, as shown in Figure 2. (Note that Abra- 
ham et al. 1996 argue that the MDS data is itself skewed to 
high asymmetry with respect to local samples.) Because the 
HDF galaxies are presumably at higher redshifts than their 
MDS counterparts, it is important to consider whether this 
shift may arise from bandshifting or fc-corrcction biases, i. e. 
whether the differences are simply a result of our observing 
ever further into the rest-frame ultraviolet for more distant 
galaxies. 

Only limited constraints are currently available for the 
redshift distribution of field galaxies to the faint limits ex- 
plored in this paper. Using arguments based on the gravi- 
tational lensing signal in well-constrained clusters. Small et 
al. (1994) and Kneib et al (1996) claim that most objects 
will lie in the range 0.5< z <2.5. In order to dotorrrnne the 
importance of bandshifting effects on galaxies at these red- 
shifts, C and A were measured for Frei et al. (1995)'s sample 
of "generic" Hubble types between T = —5 (ellipticals) and 
T — 6 (Scd galaxies) after artificially rodshifting them to 
z = 0.5, 1.0, 1.5, and 2.0. As described in Abraham et al. 
(1996), fc-corrections were applied to each image using the 
observed colour at each position on the galaxy to interpolate 
between a set of template non-evolving SEDs (kindly sup- 
plied by A. Aragon-Salamanca). The artificially rodshifted 
galaxies occupied the region of the \og{A) vs. log(C) diagram 
within the dotted polygon shown in Figure 1, which there- 
fore maps out the portion of this diagram which is expected 
to be occupied by distant non- evolving galaxies with Hubble 
T-types between —5 and 6. It appears that the majority of 
objects in the peculiar/irregular/merger sector of Figure 1 
correspond to one or more of the following: (a) to Hubble 
types later than Sd, (b) to significantly evolved intermediate 
Hubble-type spirals, or (c) to a family of objects outside the 
Hubble system (i.e. merging or peculiar galaxies). 

At this point we introduce the visual classifications 
made for the sample by RSE and vdB. These classifica- 
tions are in good agreement with classifications from the 
A — C scheme on a one-to-one basis and reproduce well the 
boundaries of the classes identified from the simulations dis- 
cussed above (Figure 1). Agreement was particularly good 
for the spheroidal systems and, where differences between 
visual and machine classifications wore found, they gener- 
ally lay at the boundaries of the 3 categories. Both vdB and 
RSE noted that beyond 7 ^ 24 mag an increasing propor- 
tion of faintest spheroidal and spiral galaxies were less reg- 
ular than brighter examples, and, in particular, there was a 
notable absence of "grand-design" spirals. Asymmetries in 
galaxies often resulted from one or more off-centred knots of 
activity, suggestive of either increased star-formation or ul- 
traviolet signal in very distant sources (and emphasizing the 
importance of accounting for bandshifting effects). A mon- 
tage of typical galaxies with 7 >23 mag is given in Figure 3 
[Plate 1]. 

Although some of the irregular galaxies in the HDF 
appear spectacularly "chain-like", the surface density of 
strictly linear systems with multiple knots and narrow trans- 
verse extents to 7 = 25 mag is ^ 2.5 arcmin"^, i.e. a 
relatively small fraction of the total number of irregulars. 
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HDF and MDS Asymmetry [oflset 0.053] 



HDF: 21.0 mag < F814W < 24.0 mag 
MDS: F814W < 22 mag 




Figure 2. Histograms of the asymmetry parameter, A, for galax- 
ies in the Hubble Deep Field (solid line) and Medium Deep Survey 
(dashed line). Note the marked skewness of the HDF data towards 
asymmetric galaxies. The MDS histogram has been normalized 
to match approximately the total counts in the HDF. The mean 
asymmetry of the HDF histogram has been reduced slightly (as 
described in the text) to correct for the effects of dithering in the 
images. 



Figure 3. (Plate 1). Montage of / >2.3 galaxies from the HDF 
with: Irr/Pcc/Mrg - columns 1 and 2; Spirals — columns 3 and 4; 
and E/SO - columns 5 and 6. Images were selected on the basis 
of the A-C classifier (see Figure 1). 



There are a larger number of objects which have a bright 

nucleus and a faint, one-sided tail, and the properties of 
these "tadpole-like" galaxies (which arc possibly related to 
the chain- like systems), will be described in a future work 
(van den Bergh et al., in preparation). 



Figure 4. (Plate 2). The U -B (top left), B -V (top right), 
V — I (bottom left), and B — I (bottom right) colours for the 
HDF sample as a function of position on the morphological log (A) 
vs. log(C) diagram. Symbols indicate each object's colour. Po- 
sitions on the morphological diagram are based upon /-band 
morphology. Sectors that subdivide approximately the "pecu- 
liar /irregular/merger", "spiral", and "carly-typc" portions of the 
diagram (defined as for Figure 1) are shown by solid lines. Note 
the strong correlation between optical colour and morphology. 
Most peculiar systems are optically blue, but many of these galax- 
ies also exhibit red U — B colours. 



of passband colours to the next is remarkably similar. The 
dominant colour axis appears horizontal, as expected from 
the close correspondence between central concentration and 
bulge-to-disk ratio (Abraham et al. 1994). The "faint blue 
galaxies" determined from optical colours correspond to low- 
central concentration, highly asymmetric objects as origi- 
nally claimed by Glazebrook et al. (1995a). Objects classified 
as early-types are predominantly red, indirectly supporting 
the accuracy of classifications made on the basis of G and 
A, and in agreement with observations suggesting that the 
early-type population in the field remains old even at faint 
flux levels. 

The distribution of irregular/peculiar/merger objects 
on the U — B diagram exhibits an interesting behaviour 
(shown more clearly in Figure 5). Nineteen of 83 objects 
in this morphological category (determined from Figure 1) 
have blue optical-near infrared colours {V — I < 0.6) and 
comparatively red UV-optical colours {U — B > —0.2). In 
fact, objects which are red inU — B tend to have high asym- 
metry irrespective of morphological class - a trend which is 
not present in other colours. They are also faint, with only 
two of the "U — B red, V — I blue" objects being brighter 
than 7 = 23 mag. One possibility is that a significant frac- 
tion of asymmetric objects have z > 2.3 corresponding to 
the Lyman discontinuity entering the U filter. Alternatively 
this may be the result of reddening from dust catalysed by 
mergers or interactions; this seems less likely as it would 
produce a trend in the other colours. 



3 COLOURS 

A very significant advantage of the HDF over earlier data 
is the extensive imaging conducted in F300W (hereafter U) , 
F450W (B) and F606W {V). In the following we will use 
instrumental magnitudes for the UBV system since our dis- 
cussion is only based on relative colours. The 7=25 mag 
limit for the present investigation not only corresponds to 
that point where sound morphological classifications can be 
determined (§2) but also to that limit at which most of the 
galaxies are reliably detected in all 4 bands; the U band is 
the shallowest exposure and is essentially complete (3 cr) at 
U ~ 27. Thus, only a few very red objects (f7 — 7 > 2) in 
our / < 25 sample arc undetected in U. 

Figure 4 [Plate 2] presents the f/-B (F300W-F450W), 
B-V (F450W - F606W), V-I (F606W - F814W), and 
B-I (F450W - F814W) colours, obtained through 0.4 arc- 
sec apertures, for the HDF sample as a function of position 
on the morphological log(^) vs. log(C) diagram. The colour 
of an individual galaxy is represented by the colour of the 
corresponding symbol on this plot. The trend from one set 



4 NUMBER COUNTS 

The high signal-to-noise ratio of the HDF data and the re- 
liability of the morphological classifications allow us to ex- 
tend the morphologically segregated galaxy counts 3 mag- 
nitudes deeper than was possible with the Medium Deep 
Survey. This dramatic improvement refiocts the fact that 
signal-to-noise is the prime limitation in image analysis; the 
increase in angular diameter distance corresponding to the 
3 magnitude improvement is probably very small. Figure 6 
combines the HDF and MDS number counts (with error 
bars based on Poisson statistics). Also shown are the no- 
evolution predictions for each type constructed as described 
in Glazebrook et al. (1995a), adopting Schecheter luminos- 
ity functions (LFs) with parameters given by Loveday et al. 
(1992), and a high normalization (<;/)* = 0.03/i^ Mpc"^). The 
elliptical panel shows the predicted counts resulting from a 
fiat faint-end slope (a = —1), instead of Loveday's original 
steep faint-end slope. The very small difference between the 
counts determined using the A — C classifications and those 
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Figure 5. The V — I vs U — B colour-colour diagram for all 
objects with / < 25 mag. Bold symbols denote objects in the 
irregular/peculiar/merger sector of Figure 1. The small number 
of objects with no U — B colour are represented by lower limits. 
Note the plume of objects with blue V — I and red U — B colours. 



determined visually by RSE and vdB arc also indicated in 
the Figure. 

The HDF counts extend the trend already identified 
by Glazebrook et al. (1995a) to much fainter magnitudes. 
The ovcr-abundancc of irrcgular/poculiar/mcrging systems 
continues to the / = 25 mag limit, at which point these 
systems represent over 40% of all galaxies. As discussed 
above, a significant fraction of these systems may lie be- 
yond a redshift z ~2.5 and there appears to be a general 
evolution of increasing asymctry with redshift. Significantly, 
some new trends emerge from Figure 6. Beyond J = 22 mag 
the spiral counts now show a significant excess over the 
no-evolution predictions. A weaker trend is seen for the 
spheroidal systems (whose counts are only marginally above 
the no-evolution prediction) and there is some evidence of a 
turn-over in the last magnitude interval. 

The interpretation of Figure 6 depends crucially on both 
the redshift distribution of the data and the faint end slope 
of the luminosity function. If the local luminosity function is 
reliably determined from recent faint surveys (c.f. Ellis et al. 
1996) and the Lyman limit is always present in distant galax- 
ies, the spiral excess would appear to imply luminosity evo- 
lution for disk galaxies in the interval 1< z <2.5. The fiat- 
tening of the spheroidal galaxy counts beyond 1 — 24 mag 
in both the visual and automated catalogues is an interest- 
ing, though currently uncertain, result. As this population 
is predominantly red, minimal evolution is implied to quite 
high redshift and steepening a would make this result yet 
tighter. The effect of curvature here is slight: in the absence 
of luminosity evolution, the bulk of the / ~24-25 spheroidals 
are expected to lie within z <2. In this respect, the decline 
in the counts is a puzzling result although it may indicate 
that the precursors of many present-day spheroidals lie in 
one of the other morphological categories. 

The large fraction of irregular/peculiar/merger systems 
detected at faint magnitudes in the HDF supports earlier 
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Figure 6. The number-magnitude relations for morphologically 
segregated samples of galaxies from the HDF and MDS. Open 
circles indicate counts obtained from automated classifications 

determined from Figure 1, closed circles indicate the results from 
the visual classifications of RSE, and crosses indicate the results 
from the visual classifications of vdB. The MDS counts arc indi- 
cated by tiie stars on each panel. The no-evolution H = 1 curves 
from Glazebrook et al. (1995a), extrapolated to 7 = 25 mag, are 
superposed. The dashed line on the E/SO diagram shows the ef- 
fect of assuming Q = 0.1. The dotted line in panel (a) shows the 
/-band number counts determined by Small et al. (1995) from 
two deep fields imaged with the Keck telescope. 



suggestions that beyond 7 ~ 24 mag the conventional Hub- 
ble system no longer provides an adequate description of the 
structural characteristics for a significant fraction of galactic 
(and possibly proto-galactic) systems. This has important 
implications for forthcoming number count and correlation 
analyses of the Hubble Deep Field. The results from fainter 
investigations may depend strongly on the prescription used 
to define galaxies amongst collections of clumpy subcompo- 
nents. Although the magnitude limit in this paper was cho- 
sen deliberately so as to minimise this ambiguity, the pecu- 
liar structure and obvious distortion of conventional galactic 
features seen in many systems at 7 < 25 mag inevitably re- 
sults in some difficulties when trying to connect these objects 
to local counterparts. At 7 < 25 mag it still seems useful 
to exploit the (approximate) connection between sectors on 
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the A — C diagram and conventional galaxy classifications 
in order to determine number counts that can be directly 
compared with those for the Hubble system. But for fainter 
galaxies it may prove more useful to model the distributions 
of C and A (or other structural parameters) directly, and 
no longer require galaxies to be directly mapped onto bins 
corresponding to local archetypes. A quantitative character- 
isation of the morphologies of faint galaxies, based on mear 
surements of physically meaningful structural parameters, 
provides an objective route forward when the conventional 
"syntEix" of morphological description, based on the Hubble 
system and references to local archetypes, has broken down. 
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